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' Abstract. Calculations of the nuclear modification factor, RdAu, for n° production 

. in dAu collisions at \/s NN = 200 GeV are presented. The applied pQCD-improved 

£f) | parton model incorporates intrinsic kj< . Nuclear multiscattering and nuclear shadowing 

are considered in the Au nucleus. Theoretical results are displayed for midrapidity and 
high pseudorapidity, 77, and compared to preliminary PHEN1X and BRAHMS data. 
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Introduction 

Midrapidity AuAu experiments at RHIC energies display strong suppression in 7r spectra 
at high p T [1. This effect was explained by jet energy loss in hot dense matter (see 
Ref. 0). To investigate whether initial state effects play a role in this suppression, 
dAu experiments were carried out. Instead of suppression, an enhancement has been 
seen in minimum bias data jSj (the "Cronin effect" in pA collisions at FERMILAB 
energies (HE]). This result indicates the absence of extra shadowing (gluon saturation), 
related to proposed wave-function modifications in the fast moving Au nucleus (see 
Ref. However, new data with centrality dependence in midrapidity dAu collisions 
[7j could provide more detailed information about the interplay between nuclear 
multiscattering [S] and (standard) nuclear shadowing 0. Analysis of recent data at 
forward rapidities jTU] can give more insight into the ^-dependence of nuclear effects. 

Here we display the results of our NLO pQCD calculations on pion production 
at 2 < pt < 10 GeV at different centralities. In parallel, we consider available data at 
forward rapidities in dAu collisions (THj, and investigate them theoretically. 

1. Calculational method in a nutshell 

We perform calculations on 7r° production in dAu collisions using a pQCD-improved 
parton model extended by a Glauber- type collision geometry [TTJ H21 EH • Introducing 
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nuclear thickness function, t Au , with a Woods -Saxon formula for Au and t^ with a 
sharp sphere for d, the invariant cross section is obtained as [14] : 

A n dAu r i r l-(X-v)/z c Hi) 

^5-= /d 2 6d 2 rt d (r)M|b-r|)l W ttt— -rr x 
d 3 p 7 s ~^ c Jvw/z c v(l-v) 

- 1 dw 



E n -%$- = / d 2 6d 2 r * d (r) t Au (\b - r|) - W 

y dz c /d 2 k Ta /d 2 k Tb / a /d(a;a,kTa,<5 2 ) h/A^b, krb, Q 2 ) x 



vw/vzc W 

3x5(1 -u>) + ^ K ab c (s, v,w,Q,Q r , Q) 

(IV 7T 
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In our next-to-leading order (NLO) calculation [TT], da/dv represents the Born cross 
section of the partonic subprocess, and K abjC (s,v,w,Q,Q r ,Q) is the corresponding 
higher order correction term ^TJ^J^I]. We fix the factorization and renormalization 
scales and connect them to the momentum of the intermediate jet, Q = Q r = (4/3)p g 
(where p q = pr/z c ), reproducing pp data at RHIC with high precision at high pr |Hj. 
The approximate form of the 3-dimensional parton distribution function (PDF) is: 

fa/ P (x a , k Ta , Q 2 ) = f a/p (x a , Q 2 ) ■ g a/p (k Ta ) . (2) 

Here, the function f a / P {x a , Q 2 ) represents the standard NLO PDF as a function of 
momentum fraction of the incoming parton, x a , at scale Q (we use MRST(cg)). 
The partonic transverse-momentum distribution is defined phenomenologically as a 2- 
dimensional Gaussian, gw^ky), and characterized by an "intrinsic h?" width fH\ ITS]. 

Nuclear multiscattering is accounted for through broadening of the incoming 
parton's /^-distribution function, namely an increase in the width of the Gaussian: 

(k 2 T ) pA = (k 2 ) pp + A(6) where A(6) = C ■ h pA (b) . (3) 

Pion production in pp collisions at RHIC energy indicates the value (fcf-)pp = 2.5 GeV 2 
|141 Ho] . Considering the multiscattering part, h p A(b) describes the number of effective 
NN collisions (or effective collision length for partons) at impact parameter b, which 
impart an average transverse momentum squared, C ^H]- At SPS energies the effectivity 
function h pA (b) was written in terms of the number of collisions suffered by the incoming 
proton in the target nucleus. For this energy region we have found a limited number of 
semihard collisions, v m — 1 = 4 and the value C = 0.35 GeV 2 [HI]; which implies a total 
broadening A ~ 1 GeV 2 . We will assume these values at RHIC energies, although the 
precise energy and rapidity dependence of v m and C is not yet verified. 

The nuclear PDF, fb/ A {xb, k^b, Q 2 ), incorporates the shadowing parametrization 
used in HIJING jH]- Extra gluon saturation (CGC) [H| would require a further 
suppression factor to be introduced. The experimental data will decide if such an extra 
suppression factor is necessary. 

Fragmentation function (NLO KKP), D^(z c ,Q 2 ), in eq. (JTJ is responsible for the 
hadronization of parton c into a pion with momentum fraction z c at scale Q = (4/3)py. 

We present results on the nuclear modification factor, defined as follows: 
1 E n dai Au /d 3 p _ E n dai Au {with nuclear effects)/d 3 j9 
dAu ~ N Un E^dal p /d?p ~ nuclear effects)/d 3 j9 ' ( ' 



Cronin Effect at Different Rapidities at RHIC 



3 



2. Results on dAu collisions at 77 = in different centrality bins 

In Fie. 1 we display the measured R^ u for tt° production at \/s NN — 200 GeV at 
different centralities: - 20%, 20 - 40%, 40 - 60% and 60 - 88% 7J. The first three 
centrality bins (up to 60 %) are similar to each other, only a slight decrease can be seen. 




Figure 1. Theoretical results on Rjau at 77 = 0, compared to PHENIX data on 7r° jTj. 



Our theoretical pQCD results on ir°, containing multiscattering and shadowing 
(solid lines), show similar tendency and overlap with the preliminary PHENIX data. 
We predicted this behavior in Ref. [Hj. These data indicate the plateau of the applied 
Woods -Saxon density profile of the Au nucleus and the applicability of a 6-independent 
shadowing in the model. The difference between data and theory in the most peripheral 
case requires further study. To indicate the importance of nuclear multiscattering in 
dAu collisions we also show results containing only shadowing (dotted lines). 

3. Minimum bias results in dAu collisions at forward rapidities 

Figure 2 summarizes BRAHMS data on charged-hadron production at forward 
pseudorapidities [12], CGC result at 77 = 3.2 in the low-pr region jH] (dashed line), 
and our pQCD results on 7r° at px > 2 GeV. We concentrate on the high-p T region, 
where no CGC result was available. The increasing strength of conventional nuclear 
shadowing with increasing pseudorapidity is indicated by pQCD results containing only 
this shadowing (dotted lines, C — 0). Including nuclear multiscattering (Cronin effect) 
we obtain reasonable agreement between data and our pQCD calculations on R^ Au (solid 
lines). Our results indicate the validity of pQCD description for dAu collisions at high 77. 
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Figure 2. Theoretical pQCD results on R^ Au for tt at 77 = 0, 1, 2.2 and 3.2, compared 
to BRAHMS data on charged particles [H3J. The CGC result is from Ref. 0. 
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